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. We derive the maximum fidelity attainable for teleportation using a shared 

£NJ , pair of d-level systems in an arbitrary pure state. This derivation provides a 

complete set of necessary and sufficient conditions for optimal teleportation 
' protocols. We also discuss the information on the teleported particle which 

is revealed in course of the protocol using a non-maximally entangled state. 
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Entanglement is a key ingredient of quantum techniques for information processing. 
One of the striking consequences of quantum entanglement is the existence of the procedure 
called quantum teleportation []]]. This procedure allows two distant parties, traditionally 
called Alice and Bob, to transmit faithfully the quantum state of a particle. The resources 
needed for this purpose is a pair of particles in a maximally entangled state shared by 
Alice and Bob, and the possibility to transmit classical messages from Alice to Bob. The 
teleportation procedure is an extremely useful tool for understanding many properties of 
quantum entanglement ||. 

An important problem in quantum information theory is the characterization of the 
entanglement exhibited by general quantum states of bipartite systems, and the evaluation 
of their capability to perform various quantum information processing tasks. In this paper 
we consider the following question. Suppose that Alice wants to teleport to Bob an unknown 
pure state of a d-leve\ particle. Alice and Bob share a single pair of <i-level particles in 
a pure state | tele) 23. What is the maximum fidelity of teleportation using such a state, and 
what teleportation protocols achieve this limit? 

We present a derivation of the maximum fidelity, which defines a set of necessary and 
sufficient conditions for a given protocol to be optimal. These conditions turn out to be 
satisfied by the standard teleportation protocol based on projections on maximally entangled 
states with appropriately adjusted Schmidt bases. The problem of optimal teleportation can 
be related to local transformations of entangled states ||: as shown by the Horodecki's, there 
is a simple algebraic link between the optimal teleportation fidelity and the maximum singlet 
fraction Q. Thus, our result provides also the maximum singlet fraction for an arbitrary 
pure state of two ci-level systems. 

Of course, use of a non-maximally entangled state makes the teleportation procedure 
imperfect. Nevertheless, we demonstrate in this paper that one can find a silver lining in such 
a case: namely, that the teleportation procedure reveals some information on the teleported 
quantum state. This information can be converted into an estimate of the quantum state 
of the particle initially possessed by Alice. We derive here an upper bound for the mean 
estimation fidelity ||, and provide an explicit recipe for constructing the quantum state 
estimate that saturates this bound. 

We shall consider a general teleportation strategy, consisting of an arbitrary measurement 
performed on Alice's side, followed by a general transformation of Bob's particle. In the 
most general case, Alice's measurement is described by a certain positive operator-valued 
measure. Such a measure can be decomposed into rank one operators, which are represented 
by projections on not necessarily normalized states |$ r )i 2 ($ r |, where the index r runs over 
all possible outcomes of Alice's measurement. The unnormalized state vector of the particle 
owned by Bob, after Alice has measured the outcome r, is given by: 

|6r)a = 12^1(1^)1® |tele) 23 ). (1) 

After having received from Alice the outcome of her measurement, Bob performs a general 
transformation of his particle, described by: 

\b r ) 3 (b r \^Y,Brs\brh(br\Bi s , (2) 

s 

where the operators B rs satisfy 
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£4E.4. = a (3) 

s 

for each r. In order to simplify the notation, we shall not write explicitly the range of the 
parameter s, which can be different for various values of r. 

We shall quantify the quality of teleportation with the help of the mean fidelity. The 
probability that Alice obtains from her measurement the outcome r is given by the scalar 
product 3(b r \b r ) 3 . The normalized state held by the Bob in this case is \b r ) 3 / \J 3 {b r \b r ) 3 . 
After the transformation of this state described by Eq. (@), its overlap with the original 
state vector \ip) is given by J2s h^BrslK)^ 2 / z{br\b r )3- Summation of this expression over 
r with the weights 3 {b r \b r ) 3) and integration over all possible input states \ip), yields the 
complete expression for the mean fidelity: 

/= /dV>£|(i 2 <$ r |® 3^1)4.(1^1® |tele) 23 )| 2 , (4) 

J rs 

where the integral / over the space of pure states is performed using the canonical measure 
invariant with respect to unitary transformations of the states \ip). 

Let us now select in the Hilbert spaces of the particles 2 and 3 the orthonormal bases 
defined by the Schmidt decomposition of the shared state | tele) 23: 

d-l 

\te\e) 23 = J2^k\k) 2 ®\k) 3 , (5) 

fc=0 

where the nonnegative real Schmidt coefficients are put in decreasing order: Ao > Ai > 
• • • Ad_i > 0. Using this basis, we may write each of the vectors |$ r )i2 as 

|^>l 2 = El0 r fc >l®|fc>2, ( 6 ) 
k=0 

where the vectors \4>r)i are not necessarily normalized. Applying this representation, the 
expression for the mean fidelity takes the form: 
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J2 x k(ip\B rs \k)((l) k r 

k=0 



(7) 



where all scalar products are now taken in a single-particle Hilbert space. This allows us 
to drop the indexes labelling the particles. The fact that the operators |$ r )i2($ r | form a 
decomposition of unity implies the following conditions on |0j?): 

£l#X#l = *«ft- (8) 

r 

The constraints imposed on the operators B rs are given by J2 S 4s4s = £-■ Our task is now 
to optimize the expression for the mean fidelity / over all possible measurements on Alice's 
side, and transformations performed by Bob. 

We shall start by deriving an upper bound on the mean fidelity of teleportation using 
the state | tele) 23. For this purpose, let us define the vectors \v%) such that: 
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EW^I = EI«JK*I- 



(9) 



fc=0 



k=0 



The vectors \u*) are uniquely defined by decomposing (4>^\ in the basis (k\ and collecting all 
the terms multiplying each of (fc|'s. The mean fidelity can be now represented as: 



E<V>|4 s k fe >W> 



fc=0 
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E E («?l^«*r.k), 

fej=0 



(10) 



where the operators Mj,- are given by the following integrals over the space of pure states 



M kl = d^immwm 



d(d + 1 



-(S kl l + \k)(l\). 



(11) 



The second explicit form of M k i is derived in the Appendix. Inserting this representation 
for the operators M k i into Eq. (]T0|), we can reduce the expression for / to the form: 
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k=0 



(12) 



The first sum over k can be transformed with the help of Eq. (|S|) multiplied by its hermitian 
conjugate: 
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EK fc k) = Tr E \k)(u k r \u l r )(l\ 

\k,l=0 , 

= Tr (E XkXil^WMi I = E^U fc >- 



vfcJ=0 



(13) 



Furthermore, with the help of the same identity given in Eq. we may convert the 
expression in the squared modulus in Eq. (12) to the form: 

J2(k\B rs \u h r ) = Tr f^-E \ u r)( k \) 

V fc=0 / 

= Tr fs P .E A fcl fc >(^ll = EA fe (0 r fe |4 

V fc=0 / fc=0 



fc=0 



Thus we have: 
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EE + E 



r k=0 
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EA fc (^|5 rs |A:) 

k=0 



(14) 



(15) 



The first term in the above expression can be easily calculated using the condition defined 
in Eq. (||), which implies that: 



E< 



Tr II = d 



and consequently the double sum over r and yields * 
the second term in Eq. ([TED we wm use the inequality 



X 2 



(16) 

d. In order to estimate 
(17) 



valid for arbitrary complex numbers Xk a - This is simply the triangle inequality for M 
complex N- dimensional vectors = (xj-x, ■ ■ ■ , XkN) with the standard quadratic norm 
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Ea=i l^fea^- When all the vectors x^ 7^ 0, the equality sign in Eq. (|I7D holds 
if and only if there exist M strictly positive numbers a\, . . . , cim such that a/x& = a^x; for 
every pair k, I. 

With the help of the triangle inequality, we can easily find an upper bound for the sum 
over rs in Eq. (p~5|): 
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fc=0 



\B rs \k) 




(18) 



The sum over rs under square root in the above expression can be estimated by 

EK^l4^>l 2 = E(' 

rs r 

<E<< 



4*>E<*l^. ft ^™l*> 



E< 



:) = d. 



(19) 



In deriving Eq. (|19|) we have implicitly assumed that 7^ 0, but of course the above 
inequalities hold also in the case when is zero. Thus we finally obtain the following 
upper bound on the mean fidelity: 



f< 



d+1 



\fe=0 / 



(20) 



We will now analyse necessary and sufficient conditions for a given teleportation protocol 
to be an optimal one. This is the case if the inequality signs in Eqs. (|18|) and fll9|) are 
replaced by equalities. Let us denote by m the maximum index for which X m is nonzero, i.e. 
A m+ i = . . . = A(i_i = 0. Of course, it is sufficient to characterize the vectors \<$) for k < m, 
and the action of the operators B rs on the subspace spanned by the vectors |0), . . . , \m). 

The inequality sign in Eq. fll8|) becomes equality if and only if there exist m + 1 nonneg- 
ative numbers ao, . . . ,a m such that: 
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ai \ k {(j> k r \Brs\k) = a*Az<#|S„IO (21) 
for any pair k,l < m. Furthermore, equality in Eq. (|T9|) takes place if and only if: 

|(0f|S„|*>r = <^|^X*|^J3„l*>- (22) 

Let us note that for a given k the scalar products (<p k \B rs \k) cannot be identically equal to 
zero. Otherwise, Eq. (E2f) implies a contradiction: 



= E M\B rs \k)\ 2 = E(M)(k\Bt s B rs \k) = = d- (23) 



Consequently, a k must be strictly positive for k < m, as discussed after Eq. (|T7|). By taking 
the squared modulus of Eq. fl2~T|), making use of Eq. (0), and performing the summation 
over s and r we find aiX k —a k Xi for each pair k, I. Thus we obtain: 

($\B r .\k) = {(f) l r \B rs \l) (24) 

for any k,l <m. Furthermore, Eq. (p2|) implies that 

A. fl |fc> =/irsfc|0") (25) 

where /z rs fc are certain complex numbers. By taking the scalar product of this identity 
with making use of Eq. (|24| ) we see that the coefficients fi rs k are independent of k: 
fJ'rsk — fJ'rs- Then by taking the scalar product of Eq. ( p5|) with the hermitian conjugated 
identity (l\Bl s = (4> l r \(i* s , and performing summation over s we obtain: 



P>ElA*r.| 2 = E<*l S J. S r.|/> = *«- (26) 
s s 

As 5Zs l/Vs| 2 — would imply that for given r all operators B rs = 0, the above identity means 
that all the vectors \(j) k ) are mutually orthogonal for k < m. Consequently, the action of the 
operator B rs in the subspace spanned by |0), . . . , \m) is equivalent, up to a multiplicative 
constant, to the action of 

m 

4 s ocEl^)W (27) 

fc=0 

Since for a given r the vectors are mutually orthogonal and have equal norm, the 
action of each of the operators B rs is proportional to the same unitary transformation on 
the relevant subspace. Thus, in order to reach the upper bound for fidelity, it is sufficient 
for Bob to perform a unitary transformation described by the following operation on the 
subspace spanned by |0), . . . \m)\ 

i m 

B r = . El#>(*l- (28) 



Necessary and sufficient conditions for Alice's measurement to be optimal are given by 
Eq. (|) and the requirement that for k < m all vectors \(j) k ) have equal norm and are 
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mutually orthogonal. It is straightforward to check that these conditions are fulfilled by the 
standard teleportation protocol Q described by 

\<P k r= P+qd ) = e 2mkp/d \(k + q)rnodd) (29) 

where p, q = 0, . . . ,d — 1 and the index r runs from to d 2 — 1. Consequently, the standard 
teleportation protocol with appropriately adjusted bases saturates the upper bound on the 
mean fidelity derived in Eq. ([20|). 

Of course, use of a nonmaximally entangled state makes the teleportation imperfect. 
However, suppose that Alice would like to use the result of her measurement to estimate the 
quantum state which has been teleported. Thus, for each outcome r of her measurement, 
she would like to assign a state {ip^ 1 ), which is her guess for the teleported state. This state 
can be represented as a result of a unitary transformation U r performed on a reference state 
|0>: 

IVD = U r \0). (30) 

Given the input state the probability that Alice's measurement yields the outcome r 
equals Y.t=o ^i\{ ( Pr\'4')\ 2 - The fidelity of the corresponding estimate is then |( , 0|'?/'^ st )| 2 = 
I^XlO)! 2 . Thus, the mean fidelity of Alice's estimate is given by: 



; cst = W d^i(vi^|o)rEA 2 fe |(0 r fe iv)l 2 (si) 

" fc=0 



Using the invariance of the measure dip with respect to unitary transformations, we may 
change the integration according to — > U r \ip). This yields: 



d-l 



/est = ££ A||<#)}| 2 |<0 r fe |[/ r 

r k=0 

= EE A 2 (tflUrMooUM), (32) 



r k=0 

where M o is defined in Eq. (p]) . By inserting its explicit form, we obtain that: 

j /d-l d-l 

\k=0 r k=0 



/est = ( E Afe + E A feE I (<f>r\Ur\0)\ 2 ) (33) 



The first double sum over k and r gives d, which follows from Eq. (|IED . The second sum can 
be estimated using the fact that for a given r all the vectors |</>Jf) are orthogonal and have 
equal norm for k < m. Thus, the second sum over k is maximized if the operator U r maps 
the vector |0) onto the subspace spanned by the vectors \4>^) corresponding to the maximum 
Afc. As Xk are ordered decreasingly, we obtain the following upper bound on the estimation 
fidelity: 

= £r (34) 
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It is straightforward to see that the optimal estimation strategy is given by 



w ' ] -^m m - (35) 



Of course, if several of have the same maximum value, then Alice can take as a guess 
any linear combination of the corresponding vectors \ip^). Let us note that the expression 
for maximum / est has analogous structure to the optimal teleportation fidelity /, with A 
replacing Efc=o -^fe- 
lt is interesting to compare two extreme cases: if the state | tele) 23 is maximally entangled, 
the maximum estimation fidelity is 1/d, which corresponds to making completely random 
guesses by Alice. This is clear, as perfect teleportation with a maximally entangled state 
cannot reveal any information on the teleported state. On the other hand, if the state | tele) 23 
is completely disentangled, the maximum estimation fidelity is 2/(d + 1), which corresponds 
to the optimal state estimation of a d- level system from a single copy |J. In this case, the 
optimal teleportation strategy reduces to the optimal state estimation procedure, with Bob 
generating on his side an imperfect copy according to the classical message obtained from 
Alice. 

In conclusion, we have derived an upper bound for fidelity of teleportation using an 
arbitrary pure bipartite system, and characterized optimal teleportation protocols. We have 
also presented an optimal strategy for estimating the quantum state given result of the 
measurement performed in course of teleportation. 

The author thanks Prof. J. H. Eberly for valuable comments on the manuscript, and 
P. Horodecki and V. Vedral for helpful discussions. This research was partially supported by 
ARO-administered MURI grant No. DAAG-19-99- 1-0125, and by NSF grant PHY-9415583. 



APPENDIX: 



In this Appendix, we evaluate integrals in Eq. (|Tl|) . Because of symmetry, it is suffi- 
cient to consider two cases: i — 0,j — 0, and i — 0,j — 1. We will use the following 
parameterization of the state vector in the basis \k): 



( e l( - cos 9 \ 
sin 9 cos (p 
Z3 sin 9 sin (p 

l Zd sin 6 sin ip 1 



(Al) 



where < £ < 2n, < 9, ip < ir/2, and z 3 , 
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_ ... _ . r ,Zd are complex numbers satisfying |z 3 | 2 + 

. + 1^1^ = 1. This parameterization is a straightforward generalization of the method used 
in 0. Following 0, the invariant volume element in this parameterization is given by: 



d^ 



(d-1) 



sin 1 



\2d-3, 



sin ip) 



2d-5 



4^-1 

xd(sin 9) d(sin (p) d£ dS 2 d~5, 



(A2) 



S 



where dS^-s is the volume element of the unit sphere S^-s- For the case % = 0,j = all 
the off-diagonals elements vanish, and we need to calculate only two elements: (0|M O o|0) = 
J dip cos 4 9 = 2/[d(d + 1)], and (l|M 00 |f) = J dip sin 2 6 cos 2 6 cos 2 tp = l/[d(d + 1)]. Due to 
symmetry, we have (A;|M o|^) = 1)] for all k ^ 0. For the operator M 01 , the only 

nonvanishing element is (0|M i|l) — J dip sin 2 9 cos 2 9 cos 2 (p = l/[d(d+ 1)]. 
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